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ABSTRACT 


Thl^Report  discusses  motion  In  the  equatorial  plane  of  an  oblate _ 
primary.  'He^-ineTude  terms  due  to  non-sphericity  through  In  the  potential. 
This  problem  Is  exactly  soluble  In  the  same  sense  that  the  ordinary  two-body 

problem  Is  exactly  soluble.  The  full  solution  Is  given,  the  connection 

T  V.,  L  5 

with  perturbation  theory  (as^^l^O)  Is  demonstrated,  and  the  stability  of 
circular  orbits  Is  discussed.  The  analytical  solution  Involves  all  three 
Incomplete  elliptic  Integrals  which  limits  Its  transparency. 
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I.  INTRODUCTION 


We  Investigate  herein  some  aspects  of  the  motion  of  a  satellite  about 
an  oblate  primary.  We  restrict  the  problem  In  two  ways.  First  we  only 
consider  motion  In  the  equatorial  plane.  Second  we  allow  departures  from 
sphericity  to  be  manifest  only  In  a  J2  term  [see  Eq.  (1)  below].  These  two 
restrictions  yield  an  analytically  tractable  problem  In  that  the  orbit  and 
the  time  dependence  of  the  location  may  be  explicitly  found.  Even  with  these 
two  restrictions  the  problem  Is  still  Important  within  the  solar  system  - 
most  natural  satellites,  be  they  moons  or  rings,  lie  In  the  equatorial  plane 
of  their  primaries  and  many  artificial  satellites  of  the  Earth  lie  In  the 
plane  of  the  celestial  equator. 

The  principal  advantage  of  being  able  to  solve  a  problem  analytically 
Is  that  one  then  may  approximate  the  solution  from  a  position  of  strength. 

For  example  expanding  sinx  In  a  power  series  In  x  can  be  rigorously  Justified, 
a  radius  of  convergence  computed,  and  with  luck  or  skill,  a  truncation  error 
established.  Here  one's  mathematical  underpinnings  are  secure.  In  contrast 
solving  for  y(x)  In 

<l^y(x)  =  -y(x)  ;  y(0)  *  0  ,  y'(0)  =  1 

Tx^ 

by  an  approximation  technique  rarely  offers  such  a  complete  set  of  luxuries. 
That  the  solution  of  this  differential  equation  Is  y(x)  »  sinx  Is  of  no  help 
unless  one  knows  that  -  precisely  when  approximation  Is  superfluous. 

We  already  know  how  to  approximately  solve  the  general  oblate  primary 
problem  via  Lagrange's  planetary  equations.  Here  we  shall  derive  some  of 
these  results  from  the  exact  solution.  Moreover,  we  shall  exploit  other 
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approximation  techniques  to  Investigate  the  orbit  and  the  time  dependence. 


It  turns  out  that  a  third  order  polynomial  In  the  distance  Is  critical 
to  determining  the  nature  of  the  solution.  When  J2  *  0  the  cubic  has  zero  for 
one  of  Its  roots  and  the  solution  for  the  motion  Is  accomplished  In  terms  of 
trigonometric  functions.  For  J2  0  the  full  cubic  must  be  considered  and 
the  solution  Involves  elliptic  functions.  Below  we  formulate  the  problem, 
develop  the  central  role  of  this  cubic,  and  analyze  Its  properties.  The 
following  sections  discuss  the  orbit,  the  time  dependence  of  the  longitude, 
and  the  time  dependence  of  the  distance.  As  appropriate  approximations  are 
used  to  rederive  the  classical  results  referred  to  above.  Finally,  because 
J2  1s  small  In  all  practical  applications,  the  notation  of  the  Keplerlan 
two-body  problem  has  been  maintained.  Of  course  for  J2  ^  0  a  Is  not  the 
semi -major  axis  nor  Is  e  the  eccentricity.  However,  as  J2  -*•  0  the  quantities 
denoted  here  by  a,  e,  etc.  do  approach  their  Keplerlan  values. 


J 
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II.  FORMULATION 

A.  Physics  and  Geometry 

We  choose  as  the  origin  the  center  of  mass  of  the  oblate  primary. 
The  equatorial  plane  of  the  primary  coincides  with  the  plane  z  -  0.  Since 
the  primary  Is  a  spheroid  there  Is  no  preferred  origin  of  longitudes. 
Cartesian  coordinates  x,  y»  and  z  are  related  to  spherical  coordinates 
r,  e,  (}>  (distance,  latitude,  and  longitude)  by 

X  -  r  cose  cos(j) 
y  *  r  cose  s1n4i 
z  =  r  sine 

We  take  this  coordinate  system  to  be  an  Inertial  one.  The  gravitational 
potential  of  the  oblate  primary,  whose  total  mass  Is  M, Is  assumed  to  be 

U(r)  *  -  ^  [1  +  (J2/2r2)(l  -3  sln^e)]  (1) 

r 

where  G  Is  the  universal  constant  of  gravitation  and  r  =  |jrJ.  The  physical 
Interpretation  of  J2  is  that  It  Is  the  difference  between  the  polar  and 
the  equatorial  moments  of  Inertia  per  unit  mass.  J2-O  herein  (If  J2  were 
negative  then  we  would  be  describing  prolate  spheroids).  Note  too  that  the 
dimensions  of  J2  are  those  of  a  length  squared. 

The  equations  of  motion  of  a  particle 
are 

r  =  -  VU(r)  (2) 


3 


IK 


where 


)»  a  (i*-rq)*cos*fi-r6*,  -re-2fb-r^*sin0cos0,  (r^2H)cos0-2r^s1n0) 


and 


J.  ^  sec0 
r  00’  r  d4»/ 


Now  dU/00  =  0  at  0  a  0  so  If  at  some  arbitrary  Instant  t  >  T  0andd 
vanish,  then  so  does  0  as  well  as  all  higher  time  derivatives  of  0.  There¬ 
fore,  with  the  initial  conditions  0(T)  =■  0(T)  =  0,  0(t)  =  0  Vt.  With  this 
simplification  the  equations  of  motion  (2)  reduce  to  (since  0U/0(t>  =  0) 

r  -  =  -aUi 


♦  2H  -  0  or  ^  .  „  (3) 

As  In  the  Keplerlan  two-body  problem  define  the  constant  L  by 

L  =  r*^  (Aa) 

and  rewrite  the  radial  equation  of  motion  as 

•r  .  LVr*  -  -  0U| 

0r 

1 0-0 

This  Is  directly  Integrable  to 

f2  .  E  -  U  -  LV(2r*)  (4b) 

I 
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where  E  is  a  constant  of  integration  easily  identified  with  the  total  energy 
per  unit  mass. 

Instead  of  dealing  with  E  and  L  we  choose  a  different  parametrization. 
Define  (p  «  6M)  a  and  e  via 


E  »  ’^/aa  .  =  paO-e^) 


(B) 


Then  since  =  -(p/r)[  1  +  J2/(2r*)]  the  conservation  of  energy  equation 
may  be  written  as 


ri  =  -yF(r) 
ar* 


(6) 


where 


F(r)  3  r*  -  2ar^  +  a*(l-e^)r  -aJg 

Note  that  for  bound  motion  E<0  so  a>0.  We  also  can  adjust  the  handedness  of 
our  coordinate  system  such  that  LiO,  therefore  l>e^.  This  is  all  we  can 
deduce,  in  general,  about  the  ranges  of  a  and  e.  However,  since  the  physi¬ 
cally  meaningful  constants  of  Integration  are  E  and  L,  our  parametrization 
in  terms  of  a,  e  is  valid  only  if  the  Jacobian  J{E,L)/(?(a,e)^0.  This  implies 
a<«-  and  ej^O. 

B.  The  Cubic  F(r) 

He  can  write  F(r)  in  a  variety  of  forms, 

F(r)  =  r*-  2ar*  +  a*(l-e*)r  -  aJ^ 

=  r  [r*-  2ar  +  a^l-e*)]  -aJg 

»  r  £r  -  a(l+e)]Cr-a(l-e)]  -aJg 
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r(r  -  r^)(r  -  r_)- 


=  (r  -  R3)(r  -  Rp)(r  -  Rq)  (7) 

where  r^  =  a(l  +  e)  and  R^^,  Rp,  Rg  are  the  three  roots  of  F.  Were  Jg  *  0 
then  the  three  roots  of  F  would  be  r^,  r_,  and  0.  Since  if  J2  vanishes 
a  and  e  are  the  seini -major  axis  and  eccentricity  it  follows  that  when  02  =  0 
r^  and  r_  are  the  maximum  and  minimum  distances  of  the  particle  from  the 
primary.  The  labeling  of  the  three  roots  of  F  when  J2  /  0  is  suggestive  of 
this  and  we  assume  R^^Rp^Rg.  We  first  show  that  F  has  three  real  roots  for 
small  J2> 

The  discriminant  of  the  cubic  2^  +a^z^  +  a22  +  ag  is  given  by  Q®+  R* 

where 

^  *  *1/3  ■  ^*2/3^^  ’  ^  “  ^*1*2  "  ^ag)/®  -  (*2/3^’ 

Here  a2  =  -2a,  a^  *  a^(l-e^),  and  ag  =  so 

q3  +  r2  =  -aVn-^  +  aNlaC-H-ge^)  +  aV  (8) 

27  27 

For  small  enough  J2(>0)  this  is  negative.  When  the  discriminant  of  a  cubic 
is  negative  then  the  cubic  has  three  real,  unequal  roots.  If  J2  increases 
from  0  until  Q*+  R*  *  0  then  the  cubic  has  three  real  roots  and  at  least  two 
of  them  are  equal.  This  corresponds  to  the  case  of  a  circular  orbit  and  will 
be  discussed  below. 

The  special  case  of  e-0,  J2  ^  0  implies  that  F  has  the  form 
F(r)  ■  r(r-a)^  -  aJ2.  In  this  case  we  would  expect  F  to  have  three  real 
roots  near  a,  a,  and  0.  Directly  from  the  analytical  solution  of  a  cubic 
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for  small  J2»  we  may  show  that  the  roots  are  approximately  equal  to  a  + 

a  -  and  0  in  this  case.  In  the  second  order  of  approximation  we  find 

(from  Newton's  method)  that  the  roots  are  approximately  equal  to  a  +  ^/5J  - 

Jg/lZa),  02/3.  The  special  case  of  J2  ~  0,  e  0  has  for  the  roots  r  +  ,  0. 

Further  specializing  to  e»0  we  find  that  the  roots  are  a,  a,  and  0. 

Returning  to  the  general  problem  we  assume  J2  is  small  enough  that 

+  R^<  0.  While  we  may  obtain  the  values  of  R.,  R_,  and  R^  directly  from 

dp  u 

the  analytical  solution  for  a  cubic  this  yields  nothing  transparent.  Instead 
let  us  compute  approximate  values  for  these  quantities  by  exploiting  the  fact 
that  as  J2-^0,  R^.  Rp,  and  Rg  must  approach  r^,  r  ,  and  0.  Using  Newton's 
method  we  find  that  the  first  order  approximations  to  the  roots  of  F  are 

Rp  -  r_  -  J2/(2er_)  <  r_ 

Rq  -  0  +  Vp  ^  ° 

where  p  =  a(l-e^).  Continuing  the  iteration  process  of  Newton's  method  we 
find,  to  second  order  in  J2  but  not  to  second  order  in  the  iteration  scheme, 

^  ^  J2/(2er+)  -  J|(l+3e)/(2er^)3 

Rp  -  r_  -  02/(2erJ  +  J2(l.3e)/(2er_)» 

Rg  0  +  J2/P  +  2J|/p  3  (10) 
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Note  that  the  J2  ^  0,  e  «  0  results  cannot  be  obtained  from  these  formulas 
merely  by  setting  e  «  0.  From  Eqs.  (9,  10)  we  conjecture  that 

Rj,  >r^>r.^Rp>Ro>0 

To  further  elucidate  the  nature  of  F(r)  we  compute  dF/dr  and  solve  for 

its  roots.  The  results  are  that  dF/dr  vanishes  if  r  =  R^;  R^  =  (a/3)  [2  + 

1/ 

(l+3e^)  '*].  From  the  values  of  d^F/dr*  at  these  two  points  we  deduce  that 
at  R^  F  has  a  local  minimum  while  at  R_  F  has  a  local  maximum.  Since 
F(+“>)=  +  “  and  F(0)  *  -aJg^O  we  further  conjecture  that 

Figure  1  contains  F(r)  for  a  =  1,  e  =  1/3,  Jg  =  1/10.  The  numerical 

values  of  Rq,  R_,  R^,  r_,R^,r^,  and  are  0.176200992,  0.281766487, 

0.398063916,  0.666666667,  1.051566846,  1.333333333,  and  1.425735091  In 

this  case.  Note  that  for  the  extrema  of  F  to  be  real  we  must  have  e*^  -1/3. 

The  results  of  this  analysis,  remembering  Eq.  (6),  Is  that  the  motion 

is  confined  to  that  range  of  r  such  that  F(r)  <  0,  i.e.  re[R_,R-].  When 

P  * 

r  =  R„  or  R_,  the  radial  velocity  vanishes  and  these  represent  the  turning 
P  O 

points  of  the  motion.  Lastly,  from  the  different  forms  of  F  In  Eq.  (7)  we 
see  that 


'^a  ^  "^p  ^  '^0  “  2a  =  +  r_ 

Ra«p  "  "  Vo  ' 


(11) 
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Having  gleaned  the  maximum  amount  of  information  concerning  F  and 
dr/dt  we  now  turn  to  a  discussion  of  the  orbit,  r(4i). 
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III.  THE  ORBIT 

A.  The  Exact  Solution 

From  the  definition  of  L  In  Eq.  (4a)  and  the  formula  for  r  In 
Eq.  (6)  we  may  derive  that  r'*  dr/d4  Is  related  to  r  =  dr/dt  via 

r  *  Lr'/r^ 

Therefore,  the  equation  of  the  orbit  Is  [p  »  a(l-e^)] 

(p.)2  =  -rFlrl  (12) 

Orient  the  xy  axes  (arbitrary  because  the  primary  Isa  spheroid) such  that  at 
t^T  r=Rp  and  ((»*0.  Suppose  that  the  earliest  time  when  r=R^  Is  t*T  +  t^ 
where  <t>  =  Let  the  next  time  when  r“R  be  t  =  T  +  t_  where  Note 

that  since  the  differential  equation  of  the  orbit  Is  even  In  «(i  we  can  build 
up  the  orbit  by  reflecting  about  (j)=0  either  one  of  these  segments.  Since 
we  know  that  Rp:£r<  R^  and  that  f*  and  (r')*  are  non>negat1ve  It  follows 
that  f  and  r'  are  >0  for  re(R  ,R  ),  t-Te(0,t.)  and  that  f  and  r'  are  <0 

K  V  a 

for  »'e(Rg,Rp),  t-Te(tg,tp).  Both  r  and  r'  vanish  If  r^R^  or  Rp.  Hence, 
the  solution  of  Eq.  (12)  Is 
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Define  the  angles  A  and  A  by 


R.-R«  r-R  R  R  -r 

a  p  u  3  P 

and  the  modulus  k(>0)  by 

K-K  Rn 

k‘  •  rr/  5®  •  kMO.1] 

''a  ''O 

Then  the  results  of  the  above  two  Integrations  are 


(13) 


(14) 


ap 


V2 

<(• 


2F(A.k) 


ap 


V2 

(^♦a) 


2F(A.k) 


where  FIs  the  Incomplete  elliptic  Integral  of  the  first  kind.  If  we 
Introduce  y(>0)  by  (note  as  J2-*-0,  k^-*-0,  1 ) 

then  these  may  be  more  compactly  written  as 


Ue  can  now  show  that  (|)p  -  2(^^.  To  see  this  observe  that  if  r  X  =  ir/2  so 
Y4»a  =  2K(k)  where  K  is  the  complete  elliptic  integral  of  the  first  kind. 

fl 

Also  we  note  that  if  r  =  R_  A  =  Tr/2  so  y(<<>„ -'!>.)  =  2K(k)  too  or  4>_  =  2(J).  =  4K(k)/Y. 

p  pa  pa 

Equation  (15)  gives  us  ^(r)  explicitly.  From  the  standard  discussions  of 
the  Keplerian  two-body  problem  we  are  used  to  dealing  with  the  inverse  of 
this,  r((^).  Since 

sn“^  (siniJj.K)  =  F(iJ;,k) 


where  sn(ij^,K)  is  the  sine  amplitude  function  of  argument  ip,  modulus  k,  we 
deduce  from  Eq.  (15),  with  u  =  y<I>/2.  that 


r(<|)) 


Rp[l-k^sn2(u,k)] 

l-msn2(u,k) 


Rpdn2(u,k) 

l-msn*(u,k) 


(16) 


where  dn(t{f,<)  is  the  delta  amplitude  function  of  argument  \p,  modulus  k  and 
m  is  an  abbreviation  for 


m  -  P 

^  “  P  _D 


=  (R  /Rg)k*,  l>m>k* 


(17) 


Equation  (16)  was  derived  from  the  upper  form  in  Eq.  (15).  Had  we  used  the 
lower  form  and  the  addition  theorem  for  sn(ijj,K),  the  identical  result  would've 
been  obtained.  This  demonstrates  the  analytical  equivalence  of  the  two 
results  and  that  Eq.  (16)  holds  over  the  entire  orbit.  It  is  clear  that  the 
orbit  is  periodic  because  the  elliptic  functions  are  periodic;  sn(ip,K)  and 
dn(4>,K)  share  4K(k)  as  a  period.  The  orbit  doesn't  repeat  however  as  can  be 
clearly  seen  in  Fig.  2.  This  is  for  a*l,  e  =  l/3,  Og'l/lO  and  shows  r(4) 
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0KB :T  KCR  64  PFRIDPS 


Fig.  2  64  orbits  for  a  =  1,  e  ■  l/3»  ■  1/10.  The  argument  of 

perlapse  advances  by  almost  180^  per  revolution  yielding  the 
double>lobed  structure. 
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for  64  complete  revolutions.  The  original  plot  was  In  three  colors,  which 
can  be  seen  In  the  gradations  of  grey.  The  values  of  k,  m,  y  and  4.  In  this 
case  are  0.603365954,  0.822443482,  0.748043275,  and  268?59733. 

B.  The  Advance  of  the  Perl apse 

A  well  known  result  from  perturbation  theory,  via  Lagrange's  equa¬ 
tions  of  motion.  Is  that  there  Is  a  secular  advance  of  the  argument  of 

perlapse.  From  Eq.  (16)  we  see  that  r(0)  ■  r(^  )  «  r(2^-)  so  that  this 

p  a 

advance  Is  evenly  distributed  over  each  half  of  the  orbit.  Its  amount  is 
A(o=  4>p-2ir»4K(k)/Y-2Tr.  To  first  order  In  Jg 

k^  -  aeJg/p*  .  y"=  l-J2(3-e)/p*  (18) 

and  as  k^-^0  K(k)  -»■  (tt/2)(1  +  kV4).  Therefore  we  calculate 

Ao)  *  37rJ2/p* 

confirming  the  perturbation  theory  value. 

C.  An  Approximate  Orbit 
As  k^  0 

sn(ij;,K)-^  s1n4^  -  (K/2)*(g;>s1n<|cosi;>)costl> 


dn(\|;,K)-4l  -  2(K/2)*s1n^tp 

Utilizing  these  approximations  and  those  In  Eq.  (18)  In  Eq.  (16)  we  are  led 
to,  after  some  algebra. 


The  first  term  is  just  the  Keplerian  orbit.  We  also  note  that,  to  first 
order  in  Jg,  r(0)  »  r(2Tr)  *  Rp  and  r(ir)»R^.  We  can  recover  the  value  of  Au 
by  remembering  that  at  an  apse  r*  *0.  Computing  this  from  Eq.  (19)  and 
assuming  (J)^  =  n +Au)/2  or  (|)p  =  2^  +  Aw,  where  Aw  is  of  order  ^2*  '"'®  iiwwdiately 
recover  the  above  value  for  Aw. 

D.  An  Alternative  Derivation 

The  approximate  form  of  r{^)  in  Eq.  (19)  can  be  derived  directly 
from  the  differential  equation  (12).  Set 

»•(<<>)  =  rQ((j))  +  r^(4.) 

where  rQ((|))  is  the  Keplerian  solution  p/(l+ecoS(t»)  and  r^  is  of  order  Jg. 
Using  this  in  Eq.  (12)  and  linearizing  in  r^  yields 

+  fr^  =  Jgg/a^ 

where 


f(<t')  =  [2rg2  -  3arQ  +  aMl-e*)]/(apr^) 

g(<(.)  =  aVQ/(2pr^) 

The  general  solution  for  r^(0)  is 


r^i^)  =  C  exp[-/^f(4)')d0’]  +  (J2/a=)exp[-/^f(^')d^']/^g(^')  • 
exp[/‘^'f(^'')d(>"]d({.’ 
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where  C  is  the  arbitrary  constant  of  integration.  After  performing  the 
integrals  we  find  that 

r, ((j))  =  Csin4i  _  J^Ccosit-f 3e  +  3e^(coS(li  ■*•  6sint)  +  e^ (Usin^<^)] 
(l+ecos<ti)2  "  2ep(Hecos4)}2 

We  set  C  =  0  because  we  don't  need  it  to  reproduce  r(0)  »r(2ir)  *R_  or  r(7r)  =R. 

p  a 

to  first  order  in  J2.  The  result  is  identical  to  that  obtained  above  from 
the  expansion  of  the  exact  solution,  Eq.  (19). 

E.  A  Processing  Ellipse 

Define  a  pseudo  semi-major  axis  and  eccentricity  a  and  e,  via 


2a  *  R.+  R„ 
a  p 

2ae  =  R.-  R_ 
a  p 

To  first  order  in  J2  they  are  given  by 
o  -  a  -  J2/(2p) 
e  »  e  +  J2(l+e^)/(2aep) 

The  function 


(20) 


^(♦1  ■  S 

represents  a  processing  ellipse  of  serai -major  axis  a  and  eccentricity  e. 

The  argument  of  peri  apse  advances  by  27r/Y  -  Z-n  which  falls  short  of  the 
actual  value  by  Trk^/2Y  «s  J2-^  0.  Hence,  although  this  form  is  intuitively 
attractive  and  correctly  indicates  the  limits  of  the  motion,  it  can't 
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adequately  replace  Eq.  (19)  as  an  approximation  to  the  orbit-even  In  lowest 
order.  In  particular  It  reproduces  all  of  Eq.  (19)  except  for  the  e*  term 
which  is  now  e*(1-<>sin<|»)  Instead  of  e*(l-s1n^(|)). 

F.  Circular  Orbits 

If  the  orbit  is  a  circle  then  r  is  a  constant.  From  Eq.  (16)  this 
implies  that  dn^(u,k)  is  proportional  to  1-msn^(u,k)  Vu.  As  dn^(<j^,ic)  » 

1  -  K^sn^(ij»,<)  this  is  possible  if  and  only  if  the  proportionality  constant 
is  unity  and  m  =  k^.  For  m  to  be  equal  to  k^  then  (i)  k^  =  0,  m  =  0  and  R,  *R„ 

a  p 

or  (ii)  k^  =  l,  m=1  and  Rp^Rg-  In  either  case  F(r)  has  a  double  root.  It 
is  a  simple  matter  to  show  from  Eq.  (6)  that  if  F(r)  =  (r  -  R^ )^(r  -  R2)  and 
one  applies  a  small  perturbation  to  the  orbit  r  =  Rp  then  for  this  orbit  to 
be  stable  it  must  be  that  R^ >  Rg.  As  unstable  orbits  are  of  no  interest  in 
this  discussion  it  follows  that  the  interesting  case  is  (i)  with  »*®Rj“Rp>0* 

From  the  analytical  solution  of  a  cubic  equation  we  know  that  a 
cubic  has  a  double  root  when  its  discriminant  vanishes.  From  Eq.  (8)  we 
see  that  forcing  Q’  +R^  to  vanish  implies  a  relationship  between  Jg.  a,  and 
e.  Since  J2  is  a  given  of  the  problem  and  circularity  is  a  statement  con¬ 
cerning  angular  momentum  rather  than  energy,  it's  really  a  relationship  for 
e(J2,a).  As  Q*  +  R^  =0  is  a  cubic  in  e*  but  a  quadratic  in  J2,  we'll  deal 
with  J2(a,  e)  here.  One  finds  that  a  vanishing  discriminant  implies 

27J2/(2aM  =  l-9e^  +  (l+Se^)*^* 

Since  J2  is  real  and  non-negative  we  must  have  e^>-1/3(cf  above).  Since  L 
is  real  and  positive  we  must  have  1  >e^  If  we  take  the  upper  sign  then  J2  is 
always  >0  for  e^e[-V3*  1]  while  if  we  take  the  lower  sign  then  J2  is  >0 
only  for  e*e[-l/3,0]. 
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If  we  actually  solve  for  the  roots  of  F(r)  -0  under  the  constraint 


Q*+  R2  =0  we  find  that  they  are  (a/3)[2±2(l +3e*)*'^*],  (a/3)[2»  (1 +3e2//2] 
(twice)  with  the  signs  the  same  order  as  above.  Our  stability  analysis 
demonstrated  that  the  physically  Interesting  solution  demanded  that  the 
double  root  be  the  larger.  This  occurs  for  the  lower  sign  and 

RQ=(2a/3)[l  -  (l+3e2)‘/='] 

Ra  »  Rp  *  (a/3)[2  +  (l+Se")*''^]  with  e2e[-l/3,0] 

Note  that  If  e^  =  -l/3  then  F  has  a  triple  root  at  2a/3  and  J2  =  8aV27. 

Figure  3  shows  F  for  e^  =  -l/9,  a  =  l.  and  32*0.107827329.  Note  also  that 
If  e^  =  0,  32*0  we  return  to  our  familiar  meaning  of  a  circular  orbit 
as  one  with  zero  eccentricity.  Remember  then  If  32*0  then  a  circular  orbit 
Is  the  one  with  maximum  angular  momentum.  Here,  with  32^0,  again  the 
angular  momentum  Is  larger  for  a  circular  orbit  (because  e^<0)  than  It  Is 
for  a  non-circular  one  (e^>0).  The  triple  root  represents  an  unstable 
orbit  but  has  the  maximum  angular  momentum. 

A  last  note:  Once  32  Is  large  enough,  Q*+R^>0  Ve*e[-1/3,1]  and  F 
has  a  single  real  root  and  a  pair  of  complex  conjugate  roots.  The  single 
real  root  represents  a  circular  orbit.  This  happens  for  As 

32  “  the  radius  of  this  unstable  orbit  (a32)*''^*  >  2a/3. 
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Fig.  3  F(r)  near  r  =  0  for  a  »  1,  »  -V9.  Jo  ’  0.107827329. 

Note  that  the  double  root  corresponds  to  a  locaf  mini mum  of  F 
and  that  this  implies  that  the  orbit  is  stable. 


IV.  THE  TIME  DEPENDENCE  OF  (J>  AND  r 

A.  The  Exact  Solution  for  tp  and  the  Period 
By  combining  Eqs.  (4a.  16)  we  have 

•j)  =  (L/Rp)[1 -msn2(u,k)]nd'*(u,k) 

where  u  =  y<P/2  and  nd(t//,K)  =  l/dn(ip,K).  The  initial  condition  is  (j)»0  when 
t  =  T.  If  we  introduce  the  incomplete  elliptic  integrals  of  the  second  and 
third  kinds,  E(ii;,K)  and  ll(ij;,n,K),  then  we  can  integrate  the  above  to  obtain 

27r(t-T)/P=  CRo-^Rp  + V^a^  ^  Rp(Ra  ' 

+  (Rp- Ro){Ra  +  Rp  +  Ro)n<e,m,k)  -  RpCR^  -  Rp)sin3cos6(l  -  k^sin^B)*/^ 

1- m  sin2  e 

[a^RpCRg-Rol^/^  (21) 

Here  P^  nAir^aVu  (e.g.  P  =  the  period  were  J2=  0).  The  angle  B  is  related  to 
via  sin6  =  sn(u,k)  or  3-am(u,k)  where  am  is  the  amplitude  function  of 

argument  t|',  modulus  k. 

Now  if  t  =  T+tg  then  4.=  (j)^  and  u  =  Y'J>a/2  =  K(k).  Note  that 
sn[K(k)k]  =  1  so  6=Tr/2.  As  the  incomplete  elliptic  integrals  F();;,ic),  E(ii^,K), 
and  n(ii^,n,K)  become  the  corresponding  complete  elliptic  integrals  K(k),  E(k), 
and  n(n,K)  when  their  argument  is  tt/2,  we  find 

Zrt/P-  lens- t 

(Rp-Ro)(R/Rp*Ro)n(i».k)|/  [a“Rp(R,-Ro)]‘''^  (22) 
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Similarly,  if  t  =  T  +  tp  then  <t>=<|>p,  u  =  Y({»p/2  *  2K(k).  But  sn[2K(k),k]  =  0 

so  b  =  7t.  As  all  of  the  incomplete  elliptic  integrals  obey  the  addition 

theorem  I(Tr,K)  =  2I(Tr/2,K),  it  follows  immediately  that  t_®2t..  Therefore, 

P  « 

not  only  is  the  orbit  periodic  but  the  time  development  of  <(>  is  periodic  too. 
Furthermore  the  period  is  tp  and  it  requires  as  much  time  for  (|)  to  increase 
from  0  to  4  as  it  does  for  (j)  to  increase  from  (j),  to  *  .  For  a  =  l,  e  =  l/3, 

o  a  p 

02  =  1/10,  2tp/p=  0.976979. 

Our  next  quest  is  the  value  of  tp  as  J2  -*-0.  From  Eq.  (22)  and  the 
fact  that  as  k^-*-0 

K(k)-^  I  (1  +< 

E(k)  7  (1  -  <  %) 

n(n.K:)-|  (l-n)-‘/2  |  1  +  (k  ^2)0  +  (1  -  n)-‘/^r^l 

one  can  show,  laboriously,  that 

tp  =  2tg  P  +  terms  of  order 

Since  =  An^aVu  it  follows  that  there  is  no  secular  change  in  a,  due  to  J2t 
to  first  order  in  J2.  This  is  a  standard  result  of  perturbation  theory. 

The  same  statement  is  true  of  the  eccentricity  and  can  be  derived  herein 
from  the  constancy  of  L  (by  definition)  and  the  absence  of  a  first  order  term 
in  P. 

B.  The  Exact  Solution  for  r 
We  can  rewrite  Eq.  (16)  as 


Remembering  the  discussion  In  ^  III  A  we  see  that 


/  C-sPrsll'/^  "  ^  t-TcCO.t^] 

\  I 

/  [-sFfs)]‘/2  =  ^ 

Ra  Tn, 

To  Integrate  this  we  need  the  following  result:  Let  P(u)  be  the  quartic 
polynomial  P  =  aQu'‘  +  a^u®  +  a2U*  +  830  +  a^,  ag^O  and  let  k  be  any  integer 
>1,  then 


2(k-1)ao 


I 


u^du  1, 


=  Zs 


k-3 


4 

+E  (2+j-2k)aj 

j  =  1 


J 


u*^'^dui. 


Here  P(s)  =  -sF(s),  ag=-l,  a^=2a,  a2  = -a*(l-e^),  a3  =  aJ2>  a^*0  and  k=2. 
Therefore  one  finds. 


=  |2[(Rp-Ro)n(X,m,k)  +  RQF(A,k)]  +  (Ra/a)[(Rp-Ro). 

Vz 

n(X,kSk)-RpF(A.k)jj/[Rp(R^-RQ)j‘/^  - 
for  t-Te[0,t,j  and 

a 

2Tr(t-T-tJ  , 

- p — ^  =  |2R3n(A,m',k)-[J2/(k^R3)J[(k"-m')F(A,k)+m’E(A.k 

/[l<p(R,-Ro)j‘*  *  * 
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for  t-Te[t,,t  J  where  tn'  = -(m-k^)/(l-ra)  =  If  one  so  desired  one 

dp  ^  P 

could  show  that  2tg^  =  tp  again  and  the  absence  of  first  order  term  in  tp/P  as 
J2  ^0.  Just  as  in  the  J2  =  0  problem  these  two  expressions  are  not  identical; 
it's  r(t),  (J)(t),  and  r(<j))  which  hold  over  both  parts  of  the  orbit  not  their 
inverses.  In  particular,  using  many  of  the  addition  theorems  for  the  elliptic 
integrals,  the  lower  form  may  be  written  as 

27r(t-T-t  ) 

- P - -  =  2[(Rp-RQ)[n(m,k)-n(X.m,k)]  +  RQ[K(k)-F(X,k)] 


+  (Ra/a)[(Rp-Ro)[^(k^k)  -^(X,k^k)]  -  Rp[K(k)  -  F(X.k)]|/ 
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